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.^ [ Abstract 

r> I We investigate the grand potential of the one-dimensional Hubbard model in 

d ' the high temperature limit, calculating the coefficients of the high temperature 

expansion (/3-expansion) of this function up to order /5^ by an alternative method. 
The results derived are analytical and do not involve any perturbation expansion 
in the hopping constant, being valid for arbitrary density of electrons in the one- 
dimensional model. In the half-filled case, we compare our analytical results for the 
specific heat and the magnetic susceptibility, in the high-temperature limit, with 
the ones obtained by Beni et al. and Takahashi's integral equations, showing that 
the latter result does not take into account the complete energy spectrum of the 
one-dimensional Hubbard model. The exact integral solution by Jiittner et al. is 
applied to the determination of the range of validity of our expansion in (3 in the 
half-filled case, for several different values of U . 

PACS numbers: 05., 05.30 -d, 05.30 Fk 

Keywords: One-Dimensional Hubbard Model, Grand Potential, High Temperature Ex- 
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1 Introduction 

The Hubbard model has been an important candidate to explain distinct physical phe- 
nomena such as itinerant magnetism in 4(i- metals |jl|, quasi one-dimensional organic salts 
(Q1D)0 and superconductivity in high T^ for two-dimensional materials^]. 

Since the earliest papers on what today is known as the Hubbard model flj, pertur- 
bation theory has been used due to the absence of exact solutions in dimensions higher 
than one. For the special case of the one-dimensional Hubbard model, Lieb and Wu@] 
applied the Bethe ansatz in order to get the analytical expression of the ground state 
wave function of the model with periodic space conditions in the half-filled case. The 
wave functions of the excited states of the Bethe ansatz and their corresponding energies 
were derived by Ovchinnikov0 from Lieb and Wu ground state wave function. That was 
the situation of the one-dimensional Hubbard model at T = in the early seventies. 

At that time, Takahashi|p derived an integral equation for the grand potential of 
the one-dimensional Hubbard model, based on the known energy spectrum of the Bethe 
ansatz solutions besides the string hypothesis. The functions that appear in such inte- 
gral equation satisfy an infinite set of coupled equations. At the same time, Shiba and 
Pincus0 numerically studied the exact thermodynamics of the Hubbard model of a one- 
dimensional model with a finite number of space sites. The longest chain included six 
sites with a periodic space boundary condition. Based on these results, they extended 
their conclusions on the behavior of the specific heat and the magnetic susceptibility to 
the thermodynamic limit. Later, Beni et al.^ applied the standard high temperature 
expansion (a perturbation theory in the hopping constant t) to derive the grand poten- 
tial of the one- dimensional Hubbard model up to order t^. The literature offers many 
examples of high temperature expansions of the grand canonical partition function for 
the Hubbard model in d— dimensions for d > 2, some of them up to order (3^ {P = -^) 



0, |Ty, |Tl|. However, all these works refer to some perturbative scheme where one of the 
characteristic constants of the model (i.e., its parameters) must be much larger than the 
other ones. 

Much more recently, the interest on theoretical aspects of the one-dimensional Hub- 
bard model returned when Shastry[^, ^, |T^ proved its integrability. A very interesting 



approach was developed where a d-dimensional quantum system at finite temperature is 
mapped onto a d + 1-dimensional classical model. In this method, the calculation of the 
grand potential of the quantum system reduces to obtaining the largest eigenvalue of the 
quantum transfer matrix of the d + 1-dimensional classical model. This was successfully 



applied to many quantum systems and, in particular, to the one-dimensional Hubbard 
model by Kliimper and Bariev in '96 for the half-filled case|T5|]. In '98, Martins and 
Ramos p6| and Jiittner et a/.|l^ fully performed the study of this model at finite tem- 
perature. In both references, the largest eigenvalue of the appropriate quantum matrix 



transfer was obtained through the Bethe ansate approach. In reference IJT^l, Jiittner et al. 
extended the results obtained by Kliimper and Bariev to any particle density. Differently 
from Takahashi's integral solution, the solution obtained by the quantum transfer matrix 
approach includes the solutions with SO (4) symmetry for the one-dimensional Hubbard 
modelpi. 



Charret et al. []19| developed a method to calculate the analytical expression of the exact 
coefficient of the high temperature expansion of the grand canonical partition function at 
each order in (3, applying this method to the calculation of the first three terms of the grand 



canonical partition function for the one-dimensional Generalized Hubbard model|]20|. The 
method is not based on the knowledge of the energy spectrum of the one-dimensional 
Hubbard model. 



In the present paper, we will apply the approach developed in reference |]T9| to get 
the coefficients associated to orders f3^ and /3^ of the high temperature expansion of the 
grand potential per site of the one-dimensional Hubbard model subject to a periodic 
space boundary condition. Our results are analytical and do not rest upon any additional 
hypothesis on the constants that characterize the model. We should point out that we do 
not perform a perturbative expansion besides the high temperature expansion, that is we 
do a /3-expansion of the grand canonical partition function for any density of electrons. 
In section 2 we present the grassmannian functions associated to the model, necessary 



to the application of the results of reference ||T9[ . In section 3 we present the coefficients 
at orders f]"^ and P^ of the /3-expansion of the grand canonical partition function of the 
one-dimensional Hubbard model. Section 4 is devoted to comparing our results to the 
ones known in the literature [p, H, 01- ^^ subsection 4.1 we compare our results to the 
perturbation expansion carried out by Beni et al. of the grand potential per site in the 
hopping term t. In subsection 4.2 we present a numerical comparison between Takahashi's 
results and ours, in the half-filled case. Finally, in subsection 4.3 we study, for different 
values of U, the range of validity of our expressions in (3, by comparing them to the 
numerical solutions provided by Jiittner et al. integral equations for the specific heat and 
magnetic susceptibility. Section 5 contains our conclusions. 



2 The One-Dimensional Hubbard Model 

The hamiltonian that describes the one-dimensional Hubbard model in the presence of 
an external constant magnetic field in the z direction is[||: 

N N 

i=\ (T=-l,l i=\ 

N 

^^bY. Y. ^aLa^,,, (1) 

i=l o-=-l,l 

where a]^ is the creation operator of an electron with spin a in the i-th site, and ajo- is the 
destruction operator of an electron with spin a in the z-th site. The first term on the r.h.s. 
of eq.(|I|) is the the hopping term of the kinetic energy operator with constant t. U is the 
strength of the interaction between electrons in the same site but with different spins. We 
have defined A^ = ^gfJ-sB, where g is the Lande's factor, //^ is the Bohr's magneton and 
B is the constant external magnetic field in the z direction. A^ is the number of space 
sites in the one-dimensional lattice. We use the convention: a =t= 1 and a =i= —1. 

The periodic boundary condition in space is implemented by imposing that ago- = ayva- 
and SLj^+i,a = aio-. Therefore, the hopping terms aj^aoo- and acj^^a^^i^^ become a[^aNa- 
and ajvo-aio- respectively. We point out that the hamiltonian ([|) is already in normal 



order, and the method by Charret et al.\19] can be applied directly. 



First of all, the high temperature expansion, j3 <^1, for the grand canonical partition 
function Z{P; /z) is 

Z{P■,^^) = Tr[e-^^] 

= Tt[1]+ Y-^ Tr[K"], (2) 

n=l '^• 

where K = H — //N, H is the hamiltonian of the system, N is the total number of 
electrons operator and // is the chemical potential with f3 = ■^; A; is the Boltzmann 
constant and T is the absolute temperature. 

We showed in reference [0] that for any self-interacting fermionic quantum system, the 
coefficients of the /^-expansion (^ can be written as multivariable Grassmann integrals. 
For the one- dimensional {d = 1), these coefficients are 

Tr[KT = / n dmdm e^'^=^ x 



x/C®(?7, r];iy = 0) IC®{f], ?]; i^ = 1) ■ ■ • IC®{f], r]-iy = n-l), 



(3) 



where rj, t] are Grassmann generators, /C® is the kernel of the K operator and the matrix 
A is given by 



A 



ATT Q 
A^T 



so that 



_^TT _ js^ii 



/^NxN —^NxN QvxAf •'■ QvxAf \ 



(4) 



(5) 



\ lATxAf QvxAT QvxAf ••• ^NxN / 

Each matrix A'^^ has dimension nA^ x nN, Inxn and flKrxAf being the identity and null 
matrices in dimension N x N, respectively. Here, A^ is the number of space sites and n 
is the power of the (3 term. The matrix A is independent of the particular model under 



consideration 19 



The kernel of the operator K for the one-dimensional Hubbard model on the lattice, 
written in terms of the Grassmann generators f]j and 77 j, is equal to 



/C® [f] 



N 

1=1 cr=±l 

N 



'^Z^Z^ HV[(l::iEln+u]N+l V[il=£ln+u]N+l+l + ^[li^n+jylTV+i V[il^n+u]N+l-li + 
1=1 a=±l 22 22 

Af 

+ / . UV(n+u)N+l Vin+u)N+l VuN+l VuN+h 
1=1 



(6) 



where the mapping p9| 

has been used. The generators ficr{xi,T„) have an equivalent mapping. These generators 
satisfy the boundary conditions: 



a) periodic boundary conditions in space: V\il:z£ln+^]N+N+i ^ ^f 



and i]r(i- 



ii^„+,,]Ar+l 



.<t) 



n+i/]N 



rj.o-al 



n+v\N+N' 



b) anti-periodic boundary condition in the temperature {u): V\(l:z:iln+n]N+i ~ ~^\ ^'^~^^ n]N v 
for / = 1, 2, ■ • ■ , A^, and a = =Fl. 

In order to write down the terms that contribute to )C®{fj, 77; u) in a simphfied way, 
we define 



N 



N 



1=1 ^ ^ 



We also define 



S{r],r];u) = J2 E{a)£{7],rj]P]a), (10) 

<T=±1 

T^{f],r];u) = J2 tT^if],v;'^;(^), (11) 

a=±l 



and 



N 
U{f],r];u) = Uj2V{n+u)N+l r]{n+u)N+l V^N+l ViyN+l, (12) 

1=1 
where E{a) = aXs — /i. Now, for the one- dimensional Hubbard model, the grassmannian 
function ]C®{f], rj; u) can be written as (see eq.(|)) 

/C®(?7, rj; v) = £{f], r]; v) + r~(?7, r]] v) + r+(?], tj] v) + W(r/, 77; v). (13) 



3 The Coefficients of the /3- Expansion of the Grand 
Potential for the One-Dimensional Hubbard Model 



In reference [Q , we calculated the coefficients of the terms of order /3^ and j3^ in eq. (^ 
for the one- dimensional Hubbard model, for arbitrary values of the constants t, [/, /i and 
B (the external magnetic field). In this section, the coefficients of orders /?^ and (3^ are 
calculated. 

The evaluation of integrals has been performed by a number of procedures (computer 
programs) developed by the authors in the symbolic system Maple 5.1. This collection 



of procedures is the computational implementation of the method described in reference 
p!9| . We have called this package ^ of procedures GINT. 

The procedure perm , contained in the package, is a useful tool to calculate the in- 
dependent non-null terms that contribute to rr[K^] and Tr[K^], and implements the 



symmetries discussed in reference pO|. The procedure gint, in its turn, calculates the 
multivariable Grassmann integrals, taking into account the property of factorization into 
sub-graphs (for details, see reference pO[|). 
We introduce a simplified notation. 



2nN 



< Ci(z/i)---C„(z/„) > = H drjjdf]! e'-'^' x 

-' 1=1 

X Oi{f],r];iyi)---Om{'n,V;^m), (14) 

which let us write the independent terms that contribute to rr[K^] in eq.(0) with n = 4 
as 

Tr[K^] =< SSSS > +4 < U££S > +2 < U£U£ > +8 < UT-T+S > + 
+ 4 < UUSS > +4 < UUUS > + < UUUU > +4 < T-£T+£ > + 
+ 8 < T-UT+£ > +4 < T-UT+U > +8 < T-T+££ > + 
+ 8 < T-T+U£ > +8 < T-T+UU > +2 < T-T+T-T+ > + 
+ 4 < T-T-T+T+ > . (15) 



In order to calculate the terms on the r.h.s. of ([T5|), we need the result of a set 
of Grassmann multivariable integrals. The procedure gint is used for obtaining such 



results. In reference [2C] we give a lengthy explanation on how to handle those integrals 



applying the results of reference [^, where multivariable Grassmann integrals are written 
as co-factors of A. 

Letting n = 5, the expression of Tr[K^] outputed by perm is 

Tr[K^] = 5 < U££££ > +5 < W^W^^ > +5 < UU£££ > + 
+5 < UU£U£ > +10 < UUT-T^£ > +5 < UUU££ > + 
+5 < UUUU£ > + < UUUUU > + < £££££ > + 
+10 < T-T+T-T+£ > +10 < r-T+T-T+U > +10 < T-UUr+£ > + 
+10 < T-T+£U£ > +10 < T-T+£££ > +10 < T-T+U££ > + 



^The package can be downloaded from the site http: /www. if. uff.hr. 



+10 < T-T+UUS > +10 < T-T+UUU > +10 < T-£T+US > + 
+10 < T-£T+£E > +10 < T-UT+EE > +10 < T'UT+UE > + 
+10 < T-Ur+UU > +10 < T+T-T-T+E > +10 < T+T-T-T+U > + 
+10 < UT-UT+E > +10 < UT-ET+E > +10 < UT'T+EE > + 
+10 < UT-T+UE > +10 < T-T-T+T+E > +10 < T-T'T+T+U > . 

(16) 

The relation between the grand potential per site W(/3; //) and the grand canonical 
partition function Z{f]; jj) is 

W(/3; /x) = - hm ^lnZ(/3;/z). (17) 



From eqs.(|l^), (0), (p^ ) and the resuhs (up to order /?^) derived in reference |^, we 
get the grand potential per site up to order /?^ for the one- dimensional Hubbard model, 

+ -e>?r}\l3^ + {-— + —UXl - -!-^2[/ + -^(72 V + 
8 "I V 64 16 " 16 16 / 

It is important to stress that the coefficients of the /3-expansion of W(/?; /i) are exact for 
any set of constants (t, f/, // and S) of this model. From eq.(^) we can obtain the strong 
limit approximation having ^ ^ t, as well as the atomic limit approximation having 
[/ <C t. No matter how large the values of the constants, the high temperature expansions 
still makes sense, provided that those values are finite. In this case, the /3— region where 
expression (|18D is bona fide is diminished. 

From expression (p!8|) , we can derive any physical quantity for the model at thermal 
equilibrium at high temperature. In the following we consider the two quantities: 

i) the specific heat at constant length and constant number of fermions Cl{(]), 



Cdf3) = -kf3-^ 







,dW{P;fi) 



dp 



(19) 



ii) the magnetic susceptibility x(/^)) 



,1 ^, d-'WiP^fi) 



xW) = -{i^g^^BY Q^,' - (20) 



In general, the available information is in terms of the density of electrons in the chain, 
instead of the chemical potential. The density of electrons is given by 

N ~ 9/i ' ^ ' 

where N^ is the number of electrons in the chain and /?, A^ and all the other constants 
of the model, are kept constants in the partial derivative. From eq. ( [I8| ) we get 



p = ^ = {--eu^ + -U lit" + — ^^u^ + -u I? - — u^ - -u^ ^^) 13^- 

^ N ^ IQ 8 ^ 384 '^ 24 ^ 768 16 ^ ^ ^ 

- {^Aue - -u \^ + - p' - - i/u + -f n- — u' + - pu^ + - ^i\^) p"" + 

^8 16 24 ^ 16 '^ A ^ 192 32 ^ ^ g P ; m ^ 

+ {-\^^U+^^U')(3' + {\p-\u)P + l. (22) 

4 Comparison to Previous Results 

In this section, we compare our results derived in section 3 to the ones presented by 
Beni et a/.[§, as well as to Takahashi's integral solution^ and to the exact integral 



solution derived by Jiitttner et al.'^j. We explore the results of these references in 
the high temperature region, where our expressions are valid, comparing two derived 
thermodynamical quantities, namely, the specific heat and the magnetic susceptibility, 
obtained by each particular approach. 

Throughout this section we have chosen t = 1; all the remaing constants in the model 
are expressed in units of t. 

4.1 Comparison to Beni et al. 

The high temperature expansion has been applied to the calculation of the coefficients 
of the /^-expansion in eq.(|D. Henderson et al. [jll| calculated the terms up to (/3t)^ 



of this series for the single-band Hubbard model in two and three-dimensional lattices. 
Bartkowiak et a/.[^ calculated the high temperature expansion of the extended Hub- 
bard model on the simple cubic lattice up to order (jSt)^. On the other hand, the high 
temperature expansion of the one-dimensional Hubbard model was carried out by Beni 
et al. [§] up to order t^ only, but to all orders in (3. We believe that it is always interesting 
to compare analytical results, mainly because in reference P] the perturbation expansion 
is done under the condition ^ <^ 1 and in the temperature region where /3t <C 1 whereas 
our results are valid for any ratio jj. 

We compare our main result (eq.(|l8D) for the grand potential per site of the one- 
dimensional Hubbard model to the expansion of eq.(8) of reference [§] up to order /?^. 
Even though both calculations allow to handle the problem with arbitrary density of 
electrons in the chain, we consider here the half-filled case (p = 1) only, when we have 
// = y. The difference between the expressions derived for the grand potential per site in 
both calculations is 

mp)-y^BPHm = ^pH\ (23) 

lb 
where yVsPHil^) is the grand potential per site derived from the expressions of Beni et 
al.. Such difference has diverse consequences for distinct physical quantities. In the case 
of the specific heats Cl{P) and Cbph{P), the difference between our result and the one 
derived by Beni et al. is 

CLif3)~CBPHiP) = l(3't\ (24) 

where Cbph{,P) is the specific heat derived from reference p. The difference between 
these two expressions is independent of [/, but the relative error decreases as U increases. 
The magnetic susceptibility is obtained from the grand potential per site W(/3) through 
eq. (pOD . We see from eq.p3|) that the difference between >V(/3) and Wbph{,P) is inde- 
pendent of the external magnetic field. Both approaches give the same expression for 
the magnetic susceptibility up to order (3^. It is worth noticing that the result for this 
thermodynamic function derived in reference [§] was obtained for ^ <C 1, whereas the 
approach presented here allows one to affirm that, up to this order in /3, it is valid for 
arbitrary values of jj. 



In reference p2[ Bartkowiak et al. asserted: " The HTSE method is exact at each order 
in the inverse temperature P" . We have a simple argument to understand this assertion. 
In the /3-expansion of the grand canonical partition function (see eq.(0)), the P"' term is 

10 



multiplied by Tr[K"], that is proportional to t'^^^^E^'^U"'^, with 2ni + n2 + ns = n. For 
even values of n, the term of highest order in the hopping constant t that contributes 
to order /3" is t". Actually, this is the only term to be calculated in order to get the 
exact coefficient of the /9-expansion of Z{j3; fi), since all the others terms proportional to 
^2ni ^^^ _ g^ __ n-j ^Q^i^ have been calculated in lower orders in the t-expansion after the 



standard high temperature expansion ||23||. When we include the term (t/?)" in the standard 
high temperature expansion, we lift the restriction ^ <C 1, under which perturbation 
theory is usually done. When n takes odd values, we do not have the term t" since we 
calculate traces in expansion (^. Therefore all the coefficients of the term /3" would have 
been calculated already in the terms proportional to t^"% (ni = 0,1,...,^^^). We only 
need to collect all the terms of the form t^"^ E"''^ [/"^ , under the condition 2ni + n2 + ns = n 
and drop the restriction ^ <^ 1. 

This analysis is coherent with result (0), where only the term of order /3^ in W(/3; /i) 
yielded a correction to the expansion obtained by Beni et a/.g. 

4.2 Comparison to Takahashi's Integral Solution 

In reference [§ , Takahashi presents a set of coupled integral equations of infinite order for 
the grand potential (see eq.(3.5a) in reference 0). We have obtained a numerical solution 
of Takahashi's equations by recursively iterating those equations up to third order. Doing 
so, we have obtained the specific heat and magnetic susceptibility in the half-filled case. 

Such iterated solution showed very good convergence: after the second or third itera- 
tions, the numerical error obtained for the grand potential was less than 1%. 

We compare our results and those coming from Takahashi's equations for two different 
values of U; namely, U = 4 and f7 = 8 (in units of t) . We have f3 ranging over the interval 
[0,0.1]. Within this range (and for the values of U under consideration) our results are 
almost exact (see section 4.3). We compare two physical quantities: the specific heat and 
the magnetic susceptibility. The specific heat curves are shown in figure |l|, for U = 4 
and U = 8. The difference between our results and Takahashi's solutions is not due to 
numerical approximations. To make this point clear, we show in figure ^ that the relative 
error between the approaches is larger than the estimated numerical error. In figure |] we 
compare the results for the magnetic susceptibility. They agree for small values of (3, but 
as f3 increases their difference again can not be explained by numerical errors. 

We conclude that, since our results can be considered exact in the given range of f3, 
something must be missing in Takahashi's result. For such small values of (3, both results 
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should agree perfectly, but that is not the case. 

4.3 Comparison to the Exact Integral Solution by Jiittner et al. 



In references [|T6| and |T^ , the thermodynamic properties of the one-dimensional Hubbard 



model are fully determined by two independent approaches; one important point is that 
solutions with built-in SO {4) symmetry appear in both, making their results exact. The 
same applies to reference ||15[ 



We take the results by Jiittner et al^Tl to discuss the validity of our expressions. 



since they have managed to write the largest eigenvalue of the suitable quantum transfer 
matrix as the solution of a few coupled integral equations. 

In reference |T^ Jiittner et al. extended the results of reference |]15| for arbitrary 



particle densities, obtaining a new set of coupled integral equations. In the half-filled 
case, even though this set is equivalent to the one presented in ||15|, it offers much better 
numerical convergence. 

We restrict the comparison of our results with the ones derived in reference Jl^ to the 



half-filled case. Our physical results are obtained by differentiating properly the grand 
potential per site in eq.(|l8D. This comparison can show the range for jS in which are 
suitable in describing the thermodynamics of the one-dimensional Hubbard model, for 
different values of U (in units of t). As in the previous subsections, we focus on the 
specific heat and the magnetic susceptibility. 

For the specific heat (see figure ^) for U = 1 a.t (3 = 0.3 the relative error between 
results is 0.63%; at /3 = 0.4 it amounts to 2,26%. For U = 4 and /5 = 0.2, the error is 
0,71%; at (3 = 0.24 it becomes 2%. For U = 8 our result at (3 = 0.1 shows an error of 
0.15% and at /3 = 0.16 it is 2%. 

For the magnetic susceptibility we get even better results, as shown in figure ^ For 
f/ = 1 and /? = 0.3 the error is 0.27%; at f3 = 0.4 the it is 0.88% (still less than 1% !). 
For ?7 = 4 and (3 = 0.3 the error is 0.39%; for (3 = 0.4 it is 1.98%. Finally, for ?7 = 8 and 
(3 = 0.2 the error is 0.3% and for f3 = 0.29 it becomes 2%. 

These comparisons stress that in all studied cases, the precision of our analytical 
solution turned out to be far much better than our initial expectations. 
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5 Conclusions 

The method developed in reference [^ can be easily applied to the one-dimensional 
Hubbard model, allowing us to derive exact analytical coefficients at each order of the 
/3-expansion of the grand potential. With the help of the procedure gint, written in 
the symbolic language Maple 5.1, the multivariable Grassmann integrals can be easily 
calculated. Besides, the property of factorization of graphs into sub-graphs, described 



in reference [Q, allows us to reduce the number of integrals to be actually calculated. 
Besides, the method does not involve any further approximation scheme. Even though 
the physics for f/ > and U < are different, the results of section 3 apply equally well 
for both cases. 

Beni et al.^ derived a perturbative expansion in the hopping constant t, valid in 
the temperature range /3t <C 1. It is simple to compare our results to theirs, since 
both methods yielded analytical results. For the sake of comparison, we consider the one- 
dimensional Hubbard model in the half-filled case. Equation ( p3D shows that the difference 
between our result for the grand potential per site and Beni et al. is proportional to t^/3'^; 
it gives a correction to the specific heat but keeps the magnetic susceptibility unchanged. 
Our result is compatible with the assertion of Bartkowiak et al. [^ that the coefficients 
derived by the standard high temperature expansion are exact at each order of /3". With 
our correction to the grand potential per site derived by Beni et al, we can drop the 
condition ^ "C 1 under which it was derived. 

Takahashi's integral solution[^] is derived from the energy spectrum of the Bethe ansatz 
solution of the one-dimensional Hubbard model plus the so-called string hypothesis. In 
the case of the specific heat, our correction to Takahashi's result is not associated with 
any numerical approximation or error. We present figure |] to show this fact. For the 
magnetic susceptibility there is a larger correction to the Takahashi's solution in the high 
temperature region. Certainly, these differences come from the fact that Takahashi's 
calculation does not take into account solutions to the one-dimensional Hubbard model 
with S0{4) symmetry p!8|]. 

In order to determine the range of validity in (3 of our analytical solution for the grand 
potential per site for different values of U (in units of t) , we considered the curves of the 
specific heat and magnetic susceptibility ioi U = 1, U = 4 and U = 8. In all cases studied 
the validity of our expressions are far beyond our initial expectations. For example, for 
U = 1 the error of our result for the magnetic susceptibility at /3 = 0.4 is less than 1%. 
The result of the standard high temperature expansion is not reliable for U = 1 and the 
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numerical solution of the coupled integral equations in the approach of Jiittner et al. [|rj] 
are very involved, while our analytical result is a very good approximation for ^ = 1 up 
to f] = 0.4. 

Finally, we should mention that the present approach opens the possibility of cal- 
culating the first terms of the /9-expansion of the grand canonical partition function of 
the Hubbard model in two space dimensions, as well as of one- dimensional models with 
inhomogeneities. We believe that improvements on the present approach will render a 
valuable tool for tackling with such problems. 
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Figure 2: Relative error in the specific heat obtained from two successive approxima- 
tions to Takahashi's equations (dashed hues) and the difference between Takahashi's and 
Charret et al. results (sohd lines) for: U = 4 and U = 8. 



Figure 3: Magnetic susceptibility from Takahashi's equations (solid line) and Charret et 
al. (dashed line) for U = 4. 



Figure 4: Specific heat from Jiittner et al. equations (solid line) and Charret et al. 
(dashed line) for: U = 1, U = 4 and U = 8. 



Figure 5: Magnetic susceptibility from Jiittner et al. equations (sohd line) and Charret 
et al. (dashed line) for: U = 1, U = 4 and U = 8. 
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